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Abstract 

We calculate the local correlation functions of charge and spin for the one-chain and 
two-chain Hubbard model using the density matrix renormalization group method and 
the recursion technique. Keeping only finite number of states we get good accuracy 
for the low energy excitations. We study the charge and spin gaps, bandwidths and 
weights of the spectra for various values of the on-site Coulomb interaction U and the 
electron filling. In the low energy part, the local correlation functions are different for 
the charge and spin. The bandwidths are proportional to t for the charge and J for 
the spin, respectively 

PACS: 71.27.+a,75.10.-b 
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1 Introduction 



The Hubbard model has played the role as a prototype of the strongly correlated electron 
system and the one-dimensional case has been well understood since the exact Bethe Ansatz 
solution by Lieb and Wu |], |], H, |, ||, ||, |7|, |§. Our main understanding of the one- dimensional 
Hubbard model is the Luttinger liquid behavior such as the charge-spin separation. The 
study for the two-dimensional Hubbard model as a higher dimensional generalization of the 
Luttinger liquid was motivated by the observation of the non-Fermi liquid behavior of the 
normal state of the high T c superconductor 0. 

In one dimension the charge and spin have different dispersion relations and gap behav- 
iors. All excitations are identified as density excitations of either charge or spin and the 
correlation functions show power-law decay. While the charge gap at half filling disappears 
as soons as we dope holes, the spin remains gapless at any filling. Some numerical calcula- 
tions for the one- dimensional Hubbard model as well as the t-J model have provided better 



understanding for the charge-spin separation in one dimension [110, [TT], O, [13], O, [la]. Jagla 



et al. observed separation of the single electron wave packet into the charge and spin density 
wave packets propagating with different velocities [fLOR . Zacher et al. identified separate 



excitations for the charge and spin in the single particle spectrum ]nj. Similarly, for the 
t-J model the dispersion energies of the charge and spin were found to scale with t and J, 
respectively |TJL |T5|j. 



The two-chain system has attracted more attention recently, related to the ladder com- 
pound [0, [17|] . The phase diagram of the two-chain Hubbard model derived by several 
authors shows diversity depending on parameters such as the on-site Coulomb interaction 



U, the electron filling n, and the interchain hopping [Tj| |TJj, [20], |21], [22[. At half filling the 
system is a spin-gapped insulator. As we dope holes the system becomes a Luttinger liquid 
for large interchain hopping as in one dimension. But for small interchain hopping the sys- 
tem becomes a spin-gapped phase with a gapless charge mode and the correlation functions 
show exponential decay. Therefore, the two-chain system shows differences from the one- 
dimensional case for small interchain hopping and at light hole doping the isotropic case, for 
which the interchain hopping is same as the intrachain hopping, belongs to this regime. The 
questions are whether there is the charge-spin separation for the two-chain system and how 
different they are from the one- dimensional case. 

The density matrix renormalization group (DMRG) method can describe the ground 
state properties and the low energy physics of a system quite accurately and it is proven 
to be more accurate in one dimension or in quasi one dimension |]2"B|, ¥2M. Hallberg studied 
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the spin dynamical correlation function of the one-dimensional Heisenberg model using the 
DMRG method and the continued fraction expansion of the Green's function, that is, the 
recursion technique p5 |. Here a certain momentum state was desired and a careful choice 
of the target states was crucial depending on the number of states kept per block to pro- 
duce accurate results for higher energy excitations and longer chains since the state does 
not remain in a given momentum sector during the DMRG iteration. However, because 
of the real space aspect of the DMRG method it is natural to consider a calculation of a 
quantity which is independent of the momentum such as the local correlation function to 
study the dynamics. In this case, one may expect better accuracy for the dynamical corre- 
lation functions with fewer states kept and target states. Later, the techniques to study the 
dynamical properties based on the DMRG method are developed by other authors. Pang 
et al. combined the DMRG method and the maximum entropy method p6| . Kiihner and 
White recently proposed an alternative method for the correction at higher frequencies [^7j . 
In this paper we study the local correlation functions of charge and spin for the one-chain 
and two-chain Hubbard model using the DMRG method and the recursion technique. We 
study the behavior of the gaps, bandwidths and weight of the spectra of charge and spin for 
various and mainly large values of U and n, then we compare the results for the one-chain 
and two- chain. 



2 Calculations 

The local correlation function such as the local density of states of electrons of a many body 



system described by the Hamiltonian f/"[28 is 



n A (i,t-t') = (0\Al(t)Mt')\0), (1) 

where |0) is the ground state of the system and Ai(t) = e lHt / h A i6 - lHt / h . j s an operator 
in coordinate space. For example, it is c\ for electron, ct-Cjj for spin, and cLcj, for charge. 
By inserting the identity J2 n \ n )( n \ — -f, where |n) are the complete set of eigenstates of H, 
the Fourier transform of Eq. (1) is 

n A {i, w) = Y,\(n\A\\0)\ 2 6(w - (E n - E )), (2) 

n 

where E is the ground state energy of the system. We define the local Green's function as 

G A (i,z) = (0\Al(z-H)- 1 A i \0). (3) 
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Then the local correlation function can be expressed as 



riA(i, w) — — ImGA(i,w + ie + E ). (4) 

7T 



The local Green's function can be calculated from the recursion technique ||29|| . In the 
Lanczos routine we choose |w ) = A^O) as the initial state. Then we get a set of orthogonal 
states which satisfy 

H\u ) = a \u ) + b^ui) (5) 

and for n > 1, 

H\u n ) = a n \u n ) + 6 n+1 |u n+ i) + 6 n |w n _i), (6) 

where a n = (u n \H\u n ) / (u n \u n ) , b = 0, and b 2 n = (u n \u n ) / ^(w n _i|u„_i) for n > 1. With the 
coefficients a's and o's above, we have a continued fraction form of Ga{x, z), 

G A (i, z) = J— . (7) 

^ - ao '-p — 

2 

2:— ai — — z _ 



In the DMRG method the system is divided into block 23, block 1 and block 4 [|23|, g4 
and the ground state can be represented as a sum of products of states in each block, 

|0>= E G iBl4l|44 |< as >|ii>|<4>. (8) 

We take the site % of Ai in the block 23 when the block 23 is in the middle of the whole 
lattice during the DMRG iteration. Then the local correlation function have little boundary 
effect for the finite size lattice although we use the open boundary condition. We get the 
initial state Aj|0) by changing only the block 23 part in |0). We prepare the ground state |0) 
from the DMRG iteration and get the coefficients a's and 6's from the recursion equations, 
Eqs. (5) and (6) until the Lanczos routine converges. 
The Hubbard model Hamiltonian^J is 

H = -tJ2 feW + H.c.) + 5>iT " TpK - \), (9) 

where t is the hopping integral and nj CT = c\ a Ci a is the number of particles with spin a at 
site i. The factor | in the second term is introduced to adjust the chemical potential for the 
particle-hole symmetry at half filling. The spin operators are 

J 3 = \ E(mr - nn), J+ = J2 4 T ^i, J- = (J + )^ (10) 
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The charge operators are 



(11) 



For the calculations of the local correlation function of spin we have 



-J 



(12) 



s,i 



which commutes with the charge operators and excites the spin only without changing the 
total charge and the local charge. For charge, if we have 



it commutes with the spin operators and excites the charge only without changing the total 
spin and the local spin. However, since it creates a up spin and a down spin at the same site, 
there is always energy cost U and we can not see the excitations in the lower band. Instead 
we have a bonding form of A Cti operator for the sites i and i + 1 (both in base 23), 



For the two-chain we choose the sites i and i + 1 in the same rung. Then we see the excitations 
in the lower band from the components of the ground state where both site i and i + 1 are 
empty. 

In order to check the accuracy of the method we calculate the local density of states of 
electrons for free case (U = 0) and compare with the exact result. In Fig. 1, we have the 
result for 1 by 20 lattice. We have the ground state as the target state during the iteration 
and the number of the states kept per block is 200. We use the open boundary condition 
for all calculations in this paper. For the DMRG result and exact result, the overall features 
of the spectra are similar and the bandwidths are the same. In particular, the low energy 
parts of the spectra are almost identical for the positions and the weights of the peaks, which 
implies great accuracy of this method for the low energy states. In this paper we calculate the 
local correlation functions of charge and spin for 1 by 32 and 2 by 16 Hubbard lattices. We 
choose parameters, U — 8, 16, 32 and n = 1.0, 0.75, 0.5. For the two-chain we have the same 
interchain hopping as the intrachain hopping. The number of states kept per block in the 
DMRG procedure is typically 200. To make the spectrum of the local correlation functions 
visible we use the Lorentzian width e = 0.05. For both charge and spin we calculate the 
counter part, the hole part with the operators A\^ and A' ^ respectively. 



■c,i 




(13) 




(14) 
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3 Results 



For the one-chain case (Figs. 2,3 and 4), at half filling the charge has a gap in the middle 
and has the particle- hole symmetry. The bandwidth of the charge is of the order of 8t. As 
we dope holes, the charge gap disappears p0[. In the particle part of the charge spectrum 
there are lower band, upper band (energy ~ U above the lower band) and another band 
(energy ~ 2U above the lower band) and this depends on weather the sites i and i + 1 are 
occupied or not before we add electrons. Therefore, the weight of the lower band reflects 
the probability that both sites i and i + 1 are empty and this weight becomes larger as the 
hole doping is larger. The border line between the particle and hole spectra is twice of the 
chemical potential (2/x) since we create two particles or two holes. As the case of small two- 
dimensional clusters |30|1, it shifts down as the hole doping is larger. Since both this quantity 
(2/i) and the gap between the lower and the upper band are of the order of U, the left edge 
of the upper band is always around 0. 

The spin does not have a gap and has the particle-hole symmetry at half filling. The 
bandwidth is of the order of 2 J (J = At 2 /U) when we take the half-width as the bandwidth. 
When we dope holes there are distinguishable inside peaks which have the same shape as 
the half filling case in a broad background. The width of the broad background is of the 
order of 8t and the width of the inside peak is proportional to J. Since there are holes which 
can move around and come back to the original position when we flip a spin, this broad 
background corresponds to the charge fluctuation and the bandwidth of this background is 
same as the bandwidth of a single electron [ 31 1 . As the hole doping is larger the weight and 
the width of the inside peak decrease. The width change approximately follows J' for the 
squeezed spin chain with hole doping, 5 = 1 — n derived by Weng et al. [§], 

J' = J[(l - 6) + sin(2vr5)/27r]. (15) 

As the on-site Coulomb interaction U increases the charge gap at half filling, which 
appears to be proportional to U, increases and the bandwidth of the charge remains same 
but the bandwidth of the spin decreases. This proves that the bandwidth of the charge scales 
with t and the bandwidth of the spin scales with J, which is consistent with the results for 



the one-dimensional t-J model fL4] , [Tq| . For large U, the distinction between the inside peak 
and outside background of the spin becomes clear. 

For the two-chain case, both charge and spin have similar features as the one-chain case. 
We show U = 32 case only in Fig. 5. Both charge and spin have the particle-hole symmetry 
at half filling. Like the one-chain case, the charge gap at half filling is proportional to U and 
disappears as we dope holes. The weight of the lower band increases as the hole doping is 
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larger. For the spin, the bandwidth at half filling is proportional to J and away from half 
filling we find the inside peaks in a broad background. The weight and the width of the 
inside peak decrease as the hole doping is larger. 

However, there are some differences from the one-chain case. The bandwidth of the 
charge at half filling is larger for the two-chain than for the one-chain because of the additional 
interchain hopping term. For the spin, the sharp edges of the inside peaks show the existence 
of the gap. The spin spectrum for n = 0.5 is significantly different from the one-chain case. 
The pseudogap feature here resembles the charge spectrum and this may be an indication of 
the absence of the charge-spin separation in this regime. 

4 Conclusions 

In this work we have studied the dynamics of charge and spin for the one-chain and two- 
chain Hubbard model. Since the DMRG method and the recursion technique produce the 
low energy part of the spectra of the local correlation functions with great accuracy, different 
behaviors of the charge and spin are clear in the low energy excitations for both one-chain 
and two-chain. The bandwidths are proportional to t for the charge and J for the spin, 
respectively. However, the background spectrum of the spin away from half filling shows 
charge behavior. The spin spectrum for the two-chain at large hole doping implies the 
different nature between one dimension and higher dimension. 
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FIGURE CAPTIONS 

Fig.l. The local density of states of electron for 1 by 20 lattice for free case (U = 0) 
with open boundary condition. The solid line corresponds to the particle part and 
the dotted line corresponds to the hole part. Top figure shows the DMRG result and 
bottom figure shows the exact result. 

Fig. 2. The local correlation functions of charge (a) and spin (b) for 1 by 32 Hubbard 
lattice with open boundary condition for U — 8. The solid line corresponds to the 
particle part and the dotted line corresponds to the hole part. These are for the half 
filling, for n = 0.75, and for n = 0.5. 

Fig.3. Same as Fig. 2 but for U = 16. 

Fig.4. Same as Fig. 2 but for U = 32. 

Fig. 5. Same as Fig. 2 but for 2 by 16 lattice and for U = 32. 
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Fig.l ( Park and Liang ) 
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Fig. 2 ( Park and Liang ) 
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Fig. 3 ( Park and Liang ) 
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Fig. 4 ( Park and Liang 
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Fig. 4 ( Park and Liang ) 
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Fig. 5 ( Park and Liang 
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Fig. 5 ( Park and Liang ) 
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